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ABSTRACT 

White dwarf (WD) binary mergers are possible progenitors to a number of unusual stars and transient phe¬ 
nomena, including type la supemovae. To date, simulations of mergers have not included magnetic fields, even 
though they are believed to play a significant role in the evolution of the merger remnant. We simulated a 0.625 

- 0.65 Mq carbon-oxygen WD binary merger in the magnetohydrodynamic moving mesh code Arepo. Each 
WD was given an initial dipole field with a surface value of ^ 10^ G. As in simulations of merging double 
neutron star binaries, we find exponential field growth within Kelvin-Helmholtz instability-generated vortices 
during the coalescence of the two stars. The final field has complex geometry, and a strength > 10^^ G at 
the center of the merger remnant. Its energy is ^ 2 x 10"^^ ergs, ^ 0.2% of the remnant’s total energy. The 
strong field likely influences further evolution of the merger remnant by providing a mechanism for angular 
momentum transfer and additional heating, potentially helping to ignite carbon fusion. 

Key words: binaries: close - magnetohydrodynamics (MHD) - stars: magnetic field - supernovae: general 

- white dwarfs 


1. INTRODUCTION 

White dwarf (WD) binaries are common end products of 
binary stellar evolution. Gravitational wave emission, mag¬ 
netic braking or the infiuence of a third body will cause a 
fraction of these to merge, producing a diversity of unusual 
stars and electromagnetic transients. In particular, for double 
carbon-oxygen (CO) WD mergers, the final outcome could 
be a m assive and rapidly rotating WD (eg. |Segretain et al. 
1997|), an accretion-induced collapse into a neutron star (NS) 


( |Saio & Nomoto|1985| ), or a nuclear explosion that might re- 


semble a type la supernova (SN la). Understanding the con¬ 
ditions that lead to each outcome requires understanding the 
merging process, which has, in the last decade, been investi¬ 
gated with increasingly sophisticated 3D hydrodynamic sim- 


ulations ([Loren- Aguilar et al.|2009[ [Pakmor et al.|2010 

Dan 

et al. 

[20121 [2014[ [Raskin et al.[[2012[ [Zhu et al.[[2013[ 

tel 

et al. 

20141. However, one fundamental piece missing in WD 


merger studies so far is magnetic fields 
Mergers (that do not immediately explode) are expected 
to produce differentially rotating merged objects - “merger 
remnants” - that are susceptible to magnetic dyna mo pro- 


Balbus 


cesses such as the magnetorotational instability (MRI; _ 

|& Hawle3^|1991| ), Tayler-Spruit dynamo (e.g. ISpruitpOn^ ), 
and the gcc dynamo (if convection occurs in the inner disk; 


_ dyr __ 

[Garcfa-Berro et al^|2Q12| ). It has therefore long been sus 
pected that they can generate stro ng fields, and re cent 2D sim¬ 
ulations of MRI in the remnant ( |Ji et al.||2013| ) have indeed 
shown amplification of a weak seed field to > 10^^ G. Mag¬ 
netic shear from these fields tran sports angular momentum 
over a timescale of ^ 10"^-10^ s ( van Kerkwijk et aL||2010 ; 
|Shen et aL]|2Q12| ) - far shorter than the thermal timescale of 
the remnant - and also (non-locally) heats the remnant. The 
latter, combined with loss of rotational support from angular 
momentum transport, could push remnant temperatures past 
the point of carbon ignition (^ 6 x 10^ K for densities be¬ 
tween 10^-10^ gcm“^), leading to either stable nuclear burn¬ 
ing or a runaway. This mechanism could potentially drive nu- 
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clear runaways even in remnants with masses below the Chan- 
drasek har Mass Mch that have t raditionally been considered 
stable dvan Kerkwijk et al.|2010| ). 

While field growth after the merger has been explored, field 
growth during the merger is also expected, and can have a pro¬ 
found impact on the post-merger magnetic evolution. Mag- 
netohydro dynamic (MHD) double NS binary merger simula¬ 
tions (eg. [Price & Rosswog]|2QQ6[ [Giacomazzo et al.||2014[ 


[Kiuchi et al.||2014| ) have found that Kelvin-Helmholtz vor¬ 

tices produced along the shear interface between the coalesc- 
ing stars can amplify field strengths by orders of magn itude 


( [Qbergaulinger et al.|2010 Zrake & MacEadyen|2013| ). T 

same should hold true for WD mergers. Motivated by this, ^ 


The 

mergers. Motivated by this, we 
present the first MHD simulation of a CO WD binary merger. 

2. METHODS 


We employ the moving-mesh code Arepo (|Springel|20101), 
which solves the equations of ideal MHD on a Voronoi mesh 
coupled with self-gravity. We operate the code in its pseudo- 
Lagrangian mode, so that the mesh-generating points that de¬ 
fine the Voronoi grid move with the local velocity of the fluid. 
To conserve angular momentum to within ^ 2% of its ini¬ 
tial value, we use the l atest improvements t o time integration 
and gradient estimate ( [Pakmor et al.||2015|). Arep o’s MHD 
implementation is desc ribed in [Pakmor et al.|(|2Ql Ij) an d [Pak 


mor & Springel[ ( [20T3 ); we use the [Powell et al.| ( [1999 ) eight- 
wave scheme for divergence control. Our simulation ignores 


Ignores 

outer hydrogen and helium layers, composition gradients, and 
nuclear reactions (negligible for sub-Mch CO WD mergers; 
Loren-Aguilar et al.|2009i[Raskin et al.|2012[ ). 

We model the merger of two CO WDs with masses of 0.625 
and 0.65 Mq, respectively, in a circular, unsynchronized bi¬ 
nary with initial separation (Tq = 2.20 x 10^ cm (corr espond¬ 
ing period Pq = 49.5 s), chosen (using the estimate of [Eggle-| 
ton[[1983[ ) such that the lower-mass WD just fills its Roche 
lobe. We chose masses typical of the narro wly peaked em¬ 
pirica l mass distribution of field CO WDs ( [Kleinman et al.[ 
2013[ ). Our initial conditions are very similar to those of the 
0.625 - 0.65 Mq binar y simulated with s moothed p article hy¬ 
drodynamics (SPH) in [Zhu et ak ( [2013 henceforth [Z1 3 [ ). As 
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Figure 1. Series of temperature T (left column) and magnetic field strength 
\B\ (right) logarithmic intensity profiles in the equatorial plane of the merger 
for four snapshots in time (rows; time indicated at the top left of each row). 

in |Z13[ both WDs are generated with a uniform initial tem¬ 
perature of 5 X 10^ K (the corresponding thermal pressure is 
dynamically irrelevant) and a uniform composition of equal 
parts carbon and oxygen by mass. They are separately re¬ 
laxed to hydrostatic equilib rium using the SPH code GASO¬ 
LINE ( [Wadsley et al.|20Q4| ) and added to Arepo by convert¬ 
ing the SPH particles to mesh-generating points while retain¬ 
ing their conservative quantities (mass, momentum and en¬ 
ergy). A uniform 10“^ gcm“^ background grid fills up a 10^^ 
cm box centered on the binary. Each WD is given a (dynam¬ 
ically irrelevant) dipole seed magnetic field with an equato¬ 
rial surface value of 10^ G (and corresponding central field of 
^ 2 X 10^ G). The fields are overlapped when the two stars 
are placed into a binary. 

The mass resolution of our simulation is mceii ~ 1 x 10~^ 
Mq . We utilize exp licit refinement and derefinement (| Vogels- 1 
[berger et al.|2012| ) to keep cell masses within a factor of two 
of 10“^ Mq and to ensure adjacent cells differ by less than a 
factor of 10 in volume. 

3. RESULTS 

We depict the evolution of the binary in Fig.f^ highlighting 
temperature T and magnetic field strength j^fT'^Fig. shows 



t{s) 

Figure 2. Total magnetic energy Eb over time, with a best fit to the rapid 
exponential growth (dashed Eb cx : Une). 

the growth of total magnetic energy Eb over time. 

In the first stage of the merger, up to ^ 180 s, the 0.625 
Mq donor WD transfers mass to the 0.65 Mq accretor for 
about 3.5 orbits before fully disrupting. Because our initial 
conditions are approximate - the WDs are not initially tidally 
deformed - mass transfer begins in spurts as the WDs stretch 
in response to the binary potential, and occurs at a rate that 
is artificially high ( [Pan et al.|2011| ). The early mass transfer 
shears the atmospheres of both WDs. As a result, Eb grows 
roughly linearly in the first ^ 100 s, reaching about quadruple 
its initial value. Since the initial mass transfer rate is artifi¬ 
cially high, this growth is likely an overestimate, but remains 
negligible compared to what follows. 

By rsj 120 s, mass transfer becomes steady, and a stream of 
material from the donor wraps around the accretor, forming 
a shear layer. Along it, the Kelvin-Helmholtz instability gen¬ 
erates a string of hot vortices that exponentially amplify their 
entrained magnetic fields. This is illustrated in Fig. [l] (row 
2), where the hot vortices along the donor-accretor interface 
correspond to regions of high field strength. At ^ 180 s, tidal 
forces between the two WDs become strong enough to fully 
disrupt the donor, which then coalesces with the accretor over 
^ 50 s. During this time, infalling donor material spirals into 
the accretor, severely deforming the accretor while carrying 
the string of magnetized vortices toward the system’s center 
of mass (CM). Fig. shows Eb growing exponentially by a 
factor of 10’ over ^ 100 s, with an ^-folding time r = 6.4 
s, comparable to the typical turnover timescale of the largest 
eddies 27r7?eddy/ Avshear ^ 3 s, where Avshear is the velocity dif¬ 
ference across the shear layer. 

By 250 s many of the vortices have merged together into 
a hot, rapidly rotating underdense void at the CM (Fig. 
row 4). Magnetic growth within the void begins to saturate as 
its magnetic and kinetic energy approach equipartition. The 
rate of field growth slows as well, with Eb growing another 
two orders of magnitude over ^ 150 s before plateauing at 
^ 2 X 10"^^ ergs at ^ 400 s. 

In Fig. [^we show the density p, T, \B\ and ratio of mag¬ 
netic to rotational energy density eBje^^x of the merger rem¬ 
nant at 400 s. The remnant consists of a dense, degeneracy- 
supported core containing ^ 60% of the remnant’s mass, a 
partly thermally supported hot envelope that surrounds the 
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Figure 3. From leftmost to rightmost column, equatorial plane (top row) and polar (bottom) logarithmic intensity profiles of density p, temperature T, magnetic 
field strength \B\ and ratio of magnetic to rotational energy density eslerot for the simulation at 400 s (~ 170 s after coalescence). The equatorial plane density 
plot includes a linear profile of the remnant core (with the same x and y scale as the logarithmic profile) to show its shape. Arrows in the magnetic field strength 
plots indicate field directions, with their lengths equal to the fraction of the field that lies along the xy plane (top frame) and xz plane (bottom). 


core, and a rotationally supported disk, a c onfigu ration simi¬ 
lar to the SPH 0.625-0.65 Mq remnant from|Z13| The Arepo 
remnant core, however, has distinctly non-axisymmetric den¬ 
sity and temperature structures, unlike the SPH simul ation 
which achieves axisymmetry ^ 170 s after coalescence ( |Z13 


online figure set Fig. 1.16). The magnetic fields are too weak 
during the merger to have an effect on merger dynamics, so 
these contrasts are due to differences in the hydrodynamic 
schemes between Arepo and SPH (C. Zhu et al. 2015, in 
preparation). 

The remnant magnetic field configuration is complex: 
while field lines are coherent along “strands” of high field 
strength, neighboring strands often point in opposite direc¬ 
tions (see Fig. [^. In the core, the volume-averaged field 
strength is 4 x G, but strands of > 10^^ G field perforate 
the core. The field remains > 10^ G near the core-disk inter¬ 
face at ~ 10^ cm, before dropping below 10^ G at > 3 x 10^ 
cm. The total magnetic field energy is ^ 0.2% the total, 
^ 0.6% the total rotational, and ^ 6% the total differential ro¬ 
tation energy of the remnant]^ This energy is roughly equally 
partitioned into toroidal and poloidal field components, with 
the ratio of poloidal to total magnetic energy Escfy/EB = 0.62. 
Studies of local field amplification within Kelvin-Helmholtz 
vortices predict magnetic growth saturates when the m agnetic 


and kinetic energy densities are close to equipar tition ( |Qber- 
Igaulinger et al.|201^|Zrake & MacFadyen|2013] ). In our sim¬ 
ulation, this is only the case for the strands of > 10^^ G field, 
where magnetic energy density is ^7% 47%) the rota¬ 

tional (differential rotation) energy density (see Fig. col¬ 
umn 4). It is possible that because the strands are distributed 
throughout the core, they drive the core’s overall evolution 
and inhibit further magnetic amplification in their surround¬ 
ings. 

Some of the magnetized accretion stream is ejected during 



Figure 4. Total magnetic energy Eb over time for the fiducial (solid blue; 
^ceii ~ 1 X 10“^ Mq, equatorial surface field strength 10^ G) simulation and 
the robustness tests. Dashed lines represent low (green; mceii ~ 5 x 10“^ 
Mq) and high resolution (red; mceii ~ 2 x 10“^ Mq) simulations. Dotted 
lines represent 1 G (magenta) and 10“^ G (cyan) low initial held simulations. 


coalescence and integrates into the inner disk ( 1-3 X 10® 

cm), producing a 10^ -10^ G field by 400 s. This field has a 
negligible hydrodynamic effect on the disk (magnetic energy 
density to pressure ratio e^/P 10“^ at 2 x 10^ cm), and, 
unlike the field in the core, has not saturated: \B\ continues to 
grow exponentially with r r\j 200 s0 

4. ROBUSTNESS TESTS 
4.1. Resolution Test 


^ Differential rotation energy of a cell is estimated with Edj-ot = 
mcc\\\v\\V xu|yjjf. 


^ The remnant disk is generally poorly resolved, even at the highest reso¬ 
lution used in Sec. |4.1[ this may artihcially slow the disk held growth rate. 
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To check that our results are not resolution-limited, we per¬ 
formed simulations, otherwise identical to the fiducial one in 
Sec. with lower and higher mass resolutions of mceii 
5 X 10“^ Mq and mceii ~ 2 x 10“^ Mq, respectively. Fig. 
compares the Eb evolution between these simulations (dashe' 
lines) and the fiducial one (solid). 

The three runs are qualitatively identical. Donor disruption 
and the start of exponential field growth occurs ^0.75 or¬ 
bital periods earlier at low resolution, and ^0.5 periods later 
at high, because the outer layers of the WDs are better cap¬ 
tured and the differences in hydrostatic equilibrium between 
Arepo and Gasoline are less pronounced at higher resolu¬ 
tion. Exponential growth rates are similar between the runs - 
the Eb ^-folding time is r = 4.9 s for the low resolution run, 
faster than r = 6.4 s for the fiducial. At high resolution, the 
growth curve appears to be separated into two phases, with 
r = 7.8 s before coalescence, and r = 3.9 s during it. The total 
magnetic energy at the end of exponential growth is also sim¬ 
ilar - at 400 s, £"5 is ^ 4 X 10^^ ergs in the low resolution run 
and ^1.5x10"^^ ergs in the high, compared to ^ 2 x 10"^^ ergs 
at the fiducial resolution. The fiducial and high resolution runs 
also qualitatively have very similar magnetic field structures 
by 400 s. Our fiducial resolution of 1 x 10“^ Mq therefore 
appears sufficient for qualitatively capturing the growth and 
final field configuration of the merger. _ 

In their MHD disk galaxy simulations, |Pakmor & Springel| 
( |2013| ) find faster field growth and higher field strength at sat- 
uration in their lowest resolution run, which they attribute to 
larger divergence errors at lower resolution. We perform a 
similar test, and also see a trend of decreasing divergence er¬ 
ror (and more accurate magnetic evolution) at higher resolu¬ 
tion, though the errors of all our simul ations are at least a 
factor of two smaller than any reported in Pakmor & Springel] 
( |2013| ). The errors are highly localized in space and trace 
steep magnetic gradients, suggesting they contribute only to 
small-scale variations in the magnetic field. 


4.2. Changing the Seed Eield Strength 

To understand the dependence of our results on the initial 
seed field, we ran two additional simulations in which we de¬ 
creased the strength of the seed field by 3 and 6 orders of 
magnitude leading to an initial equatorial surface field of 1 
G (central field ~ 2 x 10^ G), and lO'^ G (~ 20 G), respec¬ 
tively. Fig. 1^ shows their Eb evolution (dotted lines). We find 
the growth curves to be homologous between both low initial 
field runs and the fiducial one - differing only by the ratios of 
seed Eb - up to ^ 200 s, with the ^-folding time for exponen¬ 
tial amplification approximately 6.5 s for all three runs. By 
^ 250 s, the field in the fiducial simulation begins to plateau, 
while amplification (of initially weaker fields) continues for 
several hundred more seconds in the low initial field runs. For 
both runs, Eb plateaus at ^ 3 x 10"^^ ergs, comparable to the 
fiducial ~ 2 X 10"^^ ergs. Because the fields in the low initial 
field runs remain dynamically irrelevant for longer, however, 
their structures differ from that of the fiducial run and resem¬ 
ble more the crescent in Fig. row 4. The disk field does not 
saturate in any simulation - its strength is proportional to the 
strength of the seed field, and is thus much weaker in the low 
initial field runs. Our tests thus suggest that the exponential 
growth, growth timescale and plateau Eb are robust to changes 
in initial field strength, while the remnant field configuration 
is more sensitive to the choice of seed field. 


5. DISCUSSION 

We have shown that the merger of a 0.625 - 0.65 Mq CO 
WD binary produces a strong, > 10^^ G magnetic field with 
a complex structure that winds through the remnant core and 
into the inner disk. Similar to previous simulations of binary 
NS mergers, the strong field is generated by dynamo action 
within Kelvin-Helmholtz vortices formed during the coales¬ 
cence of the two WDs. Since these vortices are ubiquitous 
in WD mergers, strong magnetic fields are a likely feature of 
all merger remnants. The degree to which a field permeates 
the remnant core depends on how thoroughly the donor and 
accretor mix during coalescence, which itself is sensitive to 
initial conditions such as the degree of synchronization be¬ 
tween the WDs, or how ac curately their tidal bulg es and early 
mass transfer are captured ( |Dan et al.|2Qll|[MT^ . A parame¬ 
ter space study of magnetized mergers is needed to investigate 
the range of possible remnant field configurations. 

We note that NS mergers simulated in Eulerian grid codes 
generally show Eb growing by only a factor of 
during coalescence, compared to the ~ 10® we see, despite 
these simulations having resolutions comparable or superior 
to our low resolution Arepo run. This weaker growth is also 
inconsistent with the amplification to l ocal kinetic equipar- 
tition seen in small-scale simula tions ( [Obergaulinger et ST 
2010t|Zrake & MacFadyen|2OT^ , and is attribu ted to insum 


cient Solu t ion in the NS merg er simulatio ns ([Giacomazzo 
et al.| 2014t Kiuchi et al.||2014[ thou gh see |Dionysopoulou 
et al.|20l5] ). Giacomazzo et al.| ( |2014 ) incorp ^ted a subgria 


magnetic ampl ification model, calibrated using |Zrake & Mac 
|Fadyen| ( [2013| )s results, into their Eulerian NS merger Simula 
tion, and found Eb amplification by a factor of ^ 10^^ over a 
single dynamical time^ This suggests Arepo may be better 
able to resolve small-scale velocity structures than an Eule¬ 
rian code at comparable resolution, or better able to couple 
these structures to magnetic growth. Further work is needed 
to understand the magnetic field growth in detail. 

The density profile of the remnant remains non- 
axisymmetric for hundreds of seconds after coalescence. As 
a result, the core continues to evolve dynamically, and by 400 
s has begun to launch a pair of spiral waves into the surround¬ 
ing medium (see the density panel of Fig. [^, which trans¬ 
port angular mom entum on a timescale ri valling that for mag¬ 
netic shear. While [Kashy ap et ah] ( |20 15 1 ) report a similar spi¬ 
ral instability in their Euleria n rem nant evolution simulation, 
SPH simulations like those of |Z 13 [ rapidly achieve axisymme- 
try after coalescence and do not form spiral waves. As noted 
earlier, this difference between Arepo and SPH simulations 
is a product of the differences between their hydrodynamic 
schemes. Further study is needed to understand these differ¬ 
ences and their consequences for remnant evolution (C. Zhu 
et al. 2015, in preparation). 

The post-merger evolution of the remnant has been fol¬ 
lowed to ^ 3 X 10"^ s with axisymmet ric cylindrical (two- 
dimensional) E ulerian grid simulatio ns ([Schwab et al.||2012 


[Ji et al.[[2013 ). As described earlier, [jT et al.[ ( [2013| )s MHD 
simulation of a 0.6 - 0.6 Mq remnant shows the develop¬ 
ment of a strong magnetic field due to MRI. The subsequent 
heating and angular momentum transport due to the fields 


^ [Price & Rosswog[(|2006|s SPH simulations also show strong amplifica- 
tion; while runs using their Euler potential MHD method suffer exaggerated 
field growth fr om imprope r boundary conditions, their B-based runs show 
similar results jPrice|201^ . 
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pushes core temperatures to ignition, supporting the possibil¬ 
ity of a nuclear runaway within sub-Mch remnants. Their re¬ 
sults are, however, likely sensitive both to their initial hydro- 
dynamic conditions - which may have artificially high core 
temperatures - and their chosen seed magnetic field, a pure 
poloidal one to optimize the onset of MRL Our much stronger 
poloidal-toroidal field could substantially change post-merger 
evolution. Moreover, the persistence of a non-axisymmetric 
remnant core will lead to evolution that clearly cannot be cap¬ 
tured in a axisymmetric cylindrical grid. We therefore stress 
the need to perform high-resolution three-dimensional simu¬ 
lations of post-merger evolution to determine the final fate of 
the remnant. 

There are a number of potentially observabl e consequences 
of the magnetic fields produced by the merger. Ji et al.|p013| ) 
note the creation of a magnetized corona and biconal jet m 
their simulations, which act in concert to cause an outflow 
of material near the remnant’s poles. This outfiow eventu ally 


Beloborodov 


(2014)es 


unbinds ^ 10“'^ Mq of material, and 
timates it should lead to an optical transient with a duration of 
^ 1 day and peak L ^ 10^ Lq, which should be detectable by 
optical surveys. 

If a remnant later experiences a nuclear runaway and ex¬ 
plodes as a SN la, its magnetic field will increase the late¬ 
time emission by trapping positrons (produced by ^^Co 
decay) that would otherwise escape the ejecta. The trapping 
efficiency depends on the strength and configuration of the 
remnant magnetic field, with a locally entangled ~ 10^^ G 
field - similar t o our findings - well able to tr ap positrons 
past 1000 days ( [Ruiz-Lapu^te & Spruit|p^^ . This is in 
line with observed~late- time SN la light curves (most recently 
[Kerzendorf et al.|2014| ). 

Those remnants that do not explode will retain strong fields 
when they reach quiescence, and populate the high-mass 
tail of the distribution of isolated high-field magnetic WDs 
(|Garcia-BeTO et al.||2012| IKlilebi et ST |2013| [Wickram^ 
[inghe et al.||2014[ [Briggs et al.||2015| ). Since these remnants 
will have temperatures high enough to fuse any remaining hy¬ 
drogen and helium they possess, their properties might even- 
tually be akin to the recently discovered hot DQ WDs (e.g. 
[Dufour et al.||2013] ). These WDs have carbon-dominated 
atmospheres and Teff « 2 x 10"^ K, are often strongly mag¬ 
netized (~ 1 MG) and sometimes have monoperiodic pho¬ 
tometric variability (pos sibly due to rapid rotation). Their 
origins remain unclear (|Althaus et al.||2009t |Lawrie et al. 
201 3} Williams et al.|2013| ). Dunlap & Clemens (in press) 
recently found that, if most known hot DQs are massive 
(M > O. 95 M 0 ), their population’s velocity dispersion corre¬ 
sponds to a kinematic age much older than what would be in¬ 
ferred from their temperatures, suggesting that at least some 
hot DQs are WD-WD merger remnants. If so, their observed 
properties would constrain merger and remnant evolutionary 
models, and double-degenerate channels for SNe la. 
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